REGULARITY OF WEAK FOLIATIONS FOR THERMOSTATS 



GABRIEL P. PATERNAIN 

Abstract. Let M be a closed oriented surface endowed with a Riemannian metric 
g. We consider the flow <fi determined by the motion of a particle under the influence 
of a magnetic field and a thermostat with external field e. We show that if <fr is 
Anosov, then it has weak stable and unstable foliations of class C 1,1 if and only if 
the external field e has a global potential U, g\ := e~ 2U g has constant curvature 
and e~ u il is a constant multiple of the area form of g\. We also give necessary and 
sufficient conditions for just one of the weak foliations to be of class C 1,1 and we 
show that the combined effect of a thermostat and a magnetic field can produce an 
Anosov flow with a weak stable foliation of class C°° and a weak unstable foliation 
which is not C ' . Finally we study Anosov thermostats depending quadratically on 
the velocity and we characterize those with smooth weak foliations. In particular, 
we show that quasi-fuchsian flows as defined by Ghys in |12j can arise in this fashion. 



1. Introduction 

Let M be a closed oriented surface endowed with a Riemannian metric g and let 
SM be its unit sphere bundle. Given / 6 C°°(M, K) and e a smooth vector field on 
M we let A G C°°(SM,R) be the function given by 

X(x, v) := f(x) + (e(x),iv) 

where i indicates rotation by 7r/2 according to the orientation of the surface. The 
present paper is concerned with the dynamical system determined by the equation 

(1) -^- = A(7,7h7- 

This equation defines a flow on SM which reduces to the geodesic flow when 
e = / = 0. If A does not depend on v, i.e. e = 0, then is the magnetic flow 
associated with the magnetic field Q := f£l a , where Q a is the area form of M. The 
magnetic flow can also be seen as the Hamiltonian flow of with respect to the 

twisted symplectic form —da + tt*Q, where 7r : TM — > M is the canonical projection 
and a is the contact 1-form of the geodesic flow. If / = 0, we obtain the Gaussian 
thermostat (or isokinetic dynamics, cf. [TZj) which is reversible in the sense that 
the flip (x, v) i— > (x, —v) conjugates <p t with 0__ t (just as in the case of geodesic 
flows). Gaussian thermostats provide interesting models in nonequilibrium statistical 
mechanics El I2H] ■ 

Thus the dynamical system governed by describes the motion of a particle on M 
under the combined influence of a magnetic field fQ a and a thermostat with external 
field e. To avoid cumbersome terminology we will call simply a thermostat and we 
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will keep in mind that a magnetic force may be present. When / = 0, we will refer 
to <p as a Gaussian or pure thermostat. 

In the present paper we will study the regularity properties of the weak stable and 
unstable foliations of Anosov thermostats. Recall that the Anosov property means 
that T(SM) splits as T(SM) = E° © E u © E s in such a way that there are constants 
C > and < p < 1 < rj such that E° is spanned by the generating vector field F of 
the flow, and for all t > we have 

W-t\ E A\ < Crf- 1 and \\d<p t \ E s\\ < C p l . 

The subbundles are then invariant and Holder continuous and have smooth integral 
manifolds, the stable and unstable manifolds, which define a continuous foliation with 
smooth leaves. Also, since we are working with surfaces, the weak stable and unstable 
bundles given by 

E + = RF © E s , 
E~ = RF © E u 

are of class C 1,a , that is, they are C 1 and the transverse derivatives are a-H61der 
for some a < 1 depending on the rates of contraction and expansion of the system 
fHlE]- Moreover, if <ft is volume preserving, the work of S. Hurder and A. Katok JH] 
shows that the weak bundles have transverse derivatives which are Zygmund regular, 
and thus a-H61der for any a < 1. 

In general, thermostats are not volume preserving and in [3| it is shown that an 
Anosov thermostat on a surface preserves an absolutely continuous measure if and 
only if e has global potential, i.e., e = -VU, where U G C°°(M,R). If e = -VU, it 
is well known (cf. Section 2) that we can make a time change in <fr so that new flow 
is the magnetic flow of the metric e~ 2U g and magnetic field e~ u Q. 

The regularity of the strong bundles E s and E u for thermostats is now well un- 
derstood. Recall that the regularity of the strong bundles could change dramatically 
with time changes while the weak bundles do not change at all. Suppose that there 
exists a time change so that E s © E u is of class C 1 . Then |T%[ Theorem 2.3] im- 
plies that the reparametrized flow is a contact Anosov flow and in particular volume 
preserving. Thus e must have a global potential U and the unit sphere bundle of 
the metric e~ 2U g must be a contact type hypersurface in the symplectic manifold 
(TM, —da + Tr*e~ u Q). The question of which energy levels of a magnetic flow are 
of contact type can be decided, at least when Q is exact, in terms of Aubry-Mather 
theory. For details we refer to PjJ 122] ■ 

If we do not insist on time changes and we just ask when is E s © E u of class C 1 (or 
Lipschitz) the situation is even more rigid and it is described in [23 121 E| • The final 
result is as follows: if Q is exact, then E s © E u is of class C 1 if and only if e = Q = 
(i.e. is a geodesic flow); if Q is not exact, then E s © E u is of class C 1 if and only if 
e = 0, g has constant curvature and Q is a constant multiple of the area form. 

Thus the only question that remains to be addressed regarding regularity is that 
of the weak stable and unstable foliations. Here we show: 
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Theorem A. Let M be a closed oriented surface and let <p be an Anosov thermostat. 
Then cf) has weak stable and unstable foliations of class C 1,1 if and only if e = —VU, 
gi := e~ 2U g has constant curvature and e~ u fl is a constant multiple of the area form 

ofgi- 

Theorem A is based on results of E. Ghys [T2JE3. Theorem 4.6 in ^3] asserts that 
a smooth Anosov flow on a closed 3-manifold with weak stable and unstable foliations 
of class C 1 ' 1 , is smoothly orbit equivalent to a suspension or to what Ghys calls a 
quasi-fuchsian flow and which are described in [T2"l Theoreme B]. (In our case, since 
we are working with circles bundles the latter alternative holds.) A quasi-fuchsian 
flow ip depends on a pair of points ([<?i], [g-i\) in Teichmuller space, has smooth weak 
stable foliation C^-conjugate to the weak stable foliation of the constant curvature 
metric g\ and smooth weak unstable foliation C°°-conjugate to the weak unstable 
foliation of the constant curvature metric g 2 - Moreover, ip preserves a volume form if 
and only if [5-1] = [g 2 \. 

These results of Ghys imply that if a thermostat has weak stable and unstable 
foliations of class C 1,1 , then the Godbillon-Vey invariant of the foliations must be 
equal to 47r 2 x(M). Following Hurder and Katok [TH] and Y. Mitsumatsu [21] we 
will compute the Godbillon-Vey invariants of a thermostat and show that they are 
equal to 47r 2 x{M) if and only if the conditions in Theorem A hold. What makes the 
calculation possible is the result proved in |3] that asserts that the weak foliations 
must be transversal to the vertical fibres of SM. 

Theorem A is partially motivated by the following surprising result of M. Wo- 
jtkowski fI7\ Theorem 5.2]: if g is a metric of negative curvature on M and e a 
vector field with zero divergence (with respect to g), then the Gaussian thermostat is 
Anosov, independently of the size of e (thus the orbits of an Anosov thermostat could 
have very large geodesic curvature). In particular if g has constant negative curva- 
ture and e is the vector field dual to a harmonic 1-form, the thermostat is Anosov. 
Theorem A shows that such a flow does not have smooth weak foliations. 

Since magnetic flows are volume preserving, ^Hl Corollary 3.5] ensures that if one 
of the weak foliations is of class C 1,1 , then they are both in fact C°°. Thus Theorem 
A combined with this result implies: 

Corollary 1. Let M be a closed oriented surface and let <p be an Anosov magnetic 
flow. Then <fi has a weak foliation of class C 1 ' 1 if and only if g has constant curvature 
and Q is a constant multiple of the area form. 

Since Gaussian thermostats are reversible, a similar conclusion can be obtained for 
them in spite of the fact that they are not volume preserving: 

Corollary 2. Let M be a closed oriented surface and let be an Anosov Gaussian 
thermostat. Then has a weak foliation of class C 1,1 if and only if e = — VC7 and 
e~ 2U g has constant curvature. 

It is tempting now to speculate that in Theorem A it suffices to assume that just 
one of the weak foliations is of class C 1,1 . Surprisingly this is not the case. As we will 
see below, the combined effect of a thermostat and a magnetic field can produce an 
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Anosov flow with a weak stable foliation of class C°° and a weak unstable foliation 
which is only C 1 '" with a < 1 and no more. 

In the next result we shall assume for simplicity that e is calibrated so that its 
divergence dive = 0. As we explain in Remark 12. HI this is no restriction at all if one 
is interested only in the regularity of the weak foliations. 

Theorem B. Let M be a closed oriented surface and let <fi be an Anosov thermostat 
with dive = 0. The following are equivalent: 

(1) <p has a weak foliation of class C 1 ' 1 ; 

(2) (p has a weak foliation of class C°° ; 

(3) g has constant curvature —c 2 and there exists h G C°°(M, R) such that h 2 + 
f 2 = c 2 and 

e=^i(fVh-hVf). 

c z 

Of course, if both foliations are of class C 1,1 we must have e = V/ = in agreement 
with Theorem A. Theorem B shows that there are thermostats with a C°° weak stable 
foliation, but whose weak unstable foliation is not C 1 ' 1 . Indeed consider a closed 
hyperbolic surface (M, g) and let / be a small smooth nonconstant function on M. 
We set h := a/1 — f 2 and 

ijf 

e = 

If / is small enough the thermostat associated with (g,f,e) will be Anosov and by 
Theorem B its weak stable foliation must be of class C°°. On the other hand, its 
weak unstable foliation cannot be of class C 1 ' 1 , because otherwise by Theorem A we 
would have e = V/ = which contradicts our choice of /. 

The results above show that the dissipative quasi-fuchsian flows i/j do not appear 
as thermostats with X(x, v) = f(x) + (e(x),iv), but if we are prepared to consider 
a more general A the situation changes. Suppose q is a traceless symmetric 2-tensor 
which we also view as a function on SM. Let 

\(x, v) = f(x) + (e(x), iv) + q x (v, v) 

and consider the flow <ft on SM defined by (QJ. Again, to simplify the exposition 
we will suppose that dive = 0. Let 5 g be the divergence operator with respect to 
the metric g acting on symmetric 2-tensors. Given a symmetric 2-tensor p we can 
write p x (u,v) = (A x u,v) where A x : T X M —>■ T X M is a symmetric linear map. Set 
det g p(x) := det A x . 

Theorem C. Let M be a closed oriented surface and let <p be an Anosov thermostat 
defined by X(x,v) = f(x) + (e(x),iv) + q x (v,v), where dive = 0. Then (ft has weak 
stable and unstable foliations of class C 1 ' 1 if and only if 

e = 0, 

8 g {q-fg) = a, 

K g + det g (q + fg) = -h 2 , 
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where K g is the curvature of g and h is a non-zero constant. 

Let go be a metric with constant curvature —1 and let q be a traceless symmetric 2- 
tensor with 5 go q = 0. Such a q is the real part of a holomorphic quadratic differential, 
so they form a real vector space of dimension 6 genus(M) — 6. We will see in Sectional 
that we can find a unique metric g in the conformal class of go for which K g +det g (q) = 
—h 2 and thus if q is small enough, the thermostat with A = q is an Anosov flow with 
C°° weak foliations. We will also see that the flow is dissipative unless q = 0, so we 
obtain quasi-fuchsian flows ip which are not volume preserving. We do not know if 
all the quasi-fuchsian flows can be realized in this way, but it is quite likely that this 
is the case (see Subsection 6.2). 

2. Preliminaries 

Let M be a closed oriented surface, SM the unit sphere bundle and 7r : SM — > M 
the canonical projection. The latter is in fact a principal S^-fibration and we let V 
be the infinitesimal generator of the action of S 1 . 

Given a unit vector v G T X M, we will denote by iv the unique unit vector orthogonal 
to v such that {v, iv} is an oriented basis of T X M. There are two basic 1-forms a and 
(3 on SM which are defined by the formulas: 

<*(*,«) (0 : = (d(x,v)n(t),v); 

The form a is the canonical contact form of SM whose Reeb vector field is the 
geodesic vector field X. 

A basic theorem in 2-dimensional Riemannian geometry asserts that there exists 
a unique 1-form ip on SM (the connection form) such that i/j(V) = 1 and 

(2) da = ip A (3 

(3) d/3 — —i/j A a 

(4) dtp = -(Koix)a A (3 

where K is the Gaussian curvature of M. In fact, the form ip is given by 



Vw)(0 = (^(o),™\ 



where Z : (— e, e) — > SM is any curve with Z(0) = (x, v) and Z(0) = £ and —g- is the 
covariant derivative of Z along the curve n o Z. 

For later use it is convenient to introduce the vector field H uniquely defined by the 
conditions /3(H) = 1 and a(H) = i/)(H) = 0. The vector fields X, H and V are dual 
to a, (3 and ip and as a consequence of (j2HH) they satisfy the commutation relations 

(5) [V,X]=H, [V,H] = -X, [X,H] = KV. 

Equations (J2HIJ) also imply that the vector fields X, H and V preserve the volume 
form a A dot and hence the Liouville measure. Note that the flow of H is given by 
R 1 o g t o R, where R(x, v) = (x, iv) and g t is the geodesic flow. 
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Let A be an arbitrary smooth function on SM. For several of the results that 
we will describe below we will not need A to be a polynomial of degree < 2 in the 
velocities as in the Introduction. We may still consider a thermostat defined by (|Tj) 
and it is easy to check that 

F = X + XV 

is the generating vector field of </>. 

Now let G := —a A da = a A (3 A ip. This volume form generates the Liouville 
measure dfi of SM. 

Lemma 2.1. We have: 

(6) LpQ = V(X)&; 

(7) L H Q = 0; 

(8) L V Q = 0. 

Proof. Note that for any vector field Y, LyO = <i(iy6). Since iyO = a A (3 = 7i*fl a , 
where Q a is the area form of M, we see that Ly® = 0. Similarly, Lx® = = 0. 
Finally L F = L X Q + L XV Q = d(i xv Q) = V(X)Q. 

□ 

2.1. Time changes. Let 7(3) be a unit speed solution of 

(9) ^ = A ( 7 , 7 ')^7 / 

ds 

where X(x, v) = f(x) + (e(ar), iv) + q x (v, v). Let U G C°°(M, R) and 

t(s) := f e~ u ^dT. 

Js 

Lemma 2.2. Let 71 (t) := 7 (s(t)) and gi := e~ 2U g. Thewy^t) is a unit speed solution 
of the thermostat determined by the quadruple (gi,e u f,e 2U (e + VU),e~ u q). (Here 
V£7 is the gradient of U with respect to the metric g.) 

Proof. It is immediate to check that 71 has speed one with respect to Q\. Recall that 
the connection of g\ is given by 

D\Y = D X Y - dU(X)Y - dU(Y)X + g{X, Y)VU. 

Let us indicate derivatives with respect to t with a dot and derivatives with respect 
to s with a prime. Then we have 

° lAfl si + s 2 ^ - 2s 2 (V^/, 7 '> 7 , + s 2 VU. 



dt ds 
Using that s = e u and © we obtain: 



e 2U (fii + (e + W, ii)ij' + q(ij', 



dt 

and the lemma follows. 

□ 
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Remark 2.3. Note that we may always choose U so that div(e + VZ7) = 0. Since 
the divergence of e 2U (e + VZ7) with respect to the metric g\ is also zero, the lemma 
ensures that given a thermostat determined by (g, f, e, q) we may always perform a 
smooth time change so that the reparametrized flow is the flow of a thermostat whose 
external field has zero divergence. 

2.2. Properties of Anosov thermostats. In this subsection we collect some prop- 
erties of Anosov thermostats <fi for A arbitrary. A result of E. Ghys ^1] ensures that <fi 
is topologically conjugate to the geodesic flow of a metric of constant negative curva- 
ture and thus <fi is transitive and topologically mixing. This fact is used in the proof 
of Lemma 12.41 below. 

Recall the definition of the weak stable and unstable bundles: 

E + = RF © E s , 

E~ = RF © E u . 

Lemma 2.4. j3] For any (x,v) G SM, V(x,v) ^ E ± (x, v). 

The lemma implies that there exist unique functions r^ 1 on SM such that 

H + r + V E E+, 

H + r-V E E~. 

Note that the Anosov property implies that r + ^ r~ everywhere. Below we will need 
to use that the functions satisfy a Riccati type equation along the flow. Note that 
are as smooth as the bundles E^ (i.e. C 1,a ). 

Lemma 2.5. j3] Let r = r ± . Then 

F(r - V(X)) + r(r - V(X)) + K = 0, 
where K := K - H(X) + A 2 + F(V(X)). 

Later on we will need the following integrated version of this equation: 
Lemma 2.6. We have: 

[ {{r-V{X)f + X 2 }d^ = -^ X {M)+f [V{X)fd^. 

JSM JSM 

Proof. Since L F Q = V(A)B an easy consequence of the Stokes' theorem shows that 

/ F(r) dy, = - [ rV(X) dfi 

JSM JSM 

hence integrating the equation in Lemma f2. 51 we obtain 

-/ rV{X)dpL+ I {r-V{X))rdji = - f (K + A 2 - H(X)) dp. 

JSM JSM J SM 

Since H preserves the Liouville measure 

/ H{X) dfi = 0, 

JSM 
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and by the Gauss-Bonnet theorem 

/ Kdfi = 4n 2 X (M), 

JSM 

so the lemma follows. 

□ 

We conclude this subsection with the following simple lemma which follows right 
away from Lemma |2~31 Note that r + — r~ is function of class C l,a which never vanishes 
and without loss of generality we may assume it is always positive. 

Lemma 2.7. 

F(log(r+ - r")) = V(X) - (r+ + r~). 



3. The Godbillon-Vey invariant 

We briefly recall the definition of the Godbillon-Vey class and the Godbillon-Vey 
invariant. Let X be a smooth manifold and T a codimension-one foliation with C°° 
leaves and transversally C 2 . We suppose that the normal bundle to T is oriented so 
that there is a 1-form r whose kernel coincides with TT. By the Frobenius theorem 
there exists a 1-form rj of class C l such that dr = rj A r. The continuous 3-form 
rj A drj is closed and its cohomology class GV(T) G H 3 (X, R), called the Godbillon- 
Vey class, is an invariant of diffeomorphism and foliated concordance. When X is a 
closed oriented 3-manifold, one can also define the Godbillon- Vey invariant 

gv(T) := / rj A drj. 
Jx 

Hurder and Katok [18j proved that if T is a foliation which is transversally only 
C 1 '" for a > 1/2, then there is a natural, well-defined Godbillon-Vey invariant gv(J-), 
which extends the previous definition for JF transversally C 2 . 

The next proposition is the natural generalization to thermostats of the Mitsumatsu 
formula for geodesic flows, see [TBI Proposition 9.1] and j2"Tj . 

Proposition 3.1. Let M be a closed oriented surface and let <fi be an Anosov ther- 
mostat with A arbitrary. Let T be one of the weak foliations and suppose it is of class 
C l ' a with a > 1/2. Then 

gv{T) = 4tt 2 X (M) - 3 / {{V{\)] 2 + [V{r)] 2 ) dfi + 2 [ V{r){V 2 {\) - 2A) dfi 

JSM JSM 

where r is the unique function of class C l,a such that H + rV G TT . 

Remark 3.2. Note that if is a magnetic flow, it is volume preserving and hence 
by fHl Theorem 3.1] the weak foliations of <fi are of class C l,a for any a < 1 and 
hence gv^) is well defined. The Godbillon-Vey invariant is also well defined under 
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an appropriate pinching condition. Let 7] s be a unit vector field spanning E s and let 
r] u be a unit vector field spanning E u . Define X s and X u by the equations: 

X s (x,v,t) := hg\d(j) t (r] s (x,v))\, 

X u (x,v,t) := -log\d(f)t{r) u (x,v))\. 

We say that is r-pinched if for all t and (x, v) we have 

-X s (x,v,t) < X u (x,v,t) < rX s (x,v,t). 

T 

This condition does not depend on time changes (cf. [T2"l Section 3.3]). Note also 
that is 1-pinched if and only if <fi preserves a volume form. 

If is r-pinched with r < 2, then Corollary 1.8 in |TS] ensures that the weak 
foliations are C 1,a with a > 1/2. Thus any thermostat which is C 1 -close to a volume 
preserving flow will have a well defined Godbillon-Vey invariant. 

Proof. We shall calculate gv{F) as if T were transversally of class C 2 . This is really 
no restriction as we can always approximate the 1-form 7] by a sequence of C 1 forms 
f] n which converge to rj in the C a -topology as in the proof of Proposition 9.1 in |18j . 
In fact, the definition of the extension of the Godbillon-Vey invariant in ^H] to C 1,a 
foliations with a > 1/2 is so that the required continuity holds. 
By Lemma [2.41 V is transversal to T and hence the 1-form 

r := —Xa — r(3 + ip 

vanishes on TT and takes the value 1 on V. Hence we may take rj to be the 1-form 
given by 

r] := —iydr. 
Using the identities (j2HH) we calculate: 

77 = -( r - V(X))a + (X + V(r))(3. 

Finally we compute 

?7 a di] = {-(A + V{r)f - (r - V(X)f + (r - V{X))V{X + V{r)) 
-(A + V(r))V(r - V(X))}a A j3 A ip 

thus 

gv(F) = [ {-(A + V{r)f - (r - V{X)f + (r - V{X))V{X + V{r)) 

JSM 

-(A + V(r))V(r- V(X))}dfi. 
Using the fact that V preserves fi we can rewrite the last integral as 

gv (F) = [ {-(A + V{r)f - (r - U(A)) 2 - 2(A + V(r))V(r - V(X))} dfi. 

JSM 

Expanding and simplifying we get 

gv (T) = f {-A 2 - (r - V{X)f - 3[V(r)} 2 - AXV(r) + 2V 2 (X)(X + V(r))} d/x. 

JSM 
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Using again that V preserves \x we have 

/ \V 2 (\)dfi = - [ [U(A)] 2 d/i 

JSM JSM 



which implies 



{-A 2 - (r - V{\)) 2 - 3[V{r)} 2 + 2V{r){V 2 {\) - 2A) - 2[U(A)] 2 } dfi 



SM 



and the proposition follows from Lemma [2.61 



□ 



We will now rewrite the formula in the proposition using a bit of Fourier analysis. 
Let L 2 (SM) be the space of square integrable functions with respect to the Liouville 
measure of SM. The space L 2 (SM) decomposes into an orthogonal direct sum of 
subspaces H n , n £ Z, such that on H n , — i V is n times the identity operator (cf. 
[T4]). Given A £ C°°(SM) we can expand it as follows: 

oo 

A = $^Q fe 

fc=0 

where Qk are smooth functions such that Qk £ H-k®Hk- Set Pq := Qo and for k > 1 
set 

■{k 2 + 2). 



P k :- 



3k 2 



-V(Q k ). 



Proposition 3.3. Under the same hypothesis as in Proposition 1,9. 1\ we have 



gv{T) = Att 2 X (M) - 3 

oo 

+3£ 



k=3 



SM 

k 2 + 2^ 2 



V[r + J2 p k 



fc=0 

Q\d\i. 



dfi 



SM 



Proof. In terms of the L 2 -inner product, Proposition 13. II gives 

gv{F) = An 2 X (M) - 3 (||V(A)|| 2 + \\V(r)\\ 2 ) + 2{V(r), V 2 (\) - 2A). 
Now observe the identities 

V 2 {Q k ) = -k 2 Q k: 



-k 2 Q h 



\\V(Q k ) 



k 2 \\Q k \ 



v 2 (\)-2\ = -J2( 2 + k2 )Qk- 



k=0 

Using that elements in H_ k © H k are orthogonal to elements in © Hi for k ^ I it 
is now straightforward to check that the identity claimed in the proposition holds. 

□ 

We now specialize Proposition 13.31 to the case that we are interested in, i.e., when 
A = / + V{6) + q, where 9 is the 1-form dual to the external field E and q £ H^2 © H2- 
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Corollary 3.4. Let M be a closed oriented surface and let be an Anosov thermostat 
defined by A = / + V(9) + q. Let T be one of the weak foliations and suppose it is of 
class C 1,a with a > 1/2. Then 

gv{T) = 4tt 2 x (M) - 3 / [V{r + 9 - V(q)/2)] 2 dfi 

J SM 

where r is the unique function of class C 1,Q such that H + rV G TJF. 
Proof. It suffices to note that since Qi = V{9) and Qi = q, then 

Pi = -V{V{9)) = 9 

and 

P 2 = -V(q)/2. 

□ 

4. Proof of Theorem A 

Suppose first that e = — VZ7, g± := e~ 2U g has constant curvature and e~ u Q is a 
constant multiple of the area form of g\. On account of Lemma 12.21 we can make a 
smooth time change to the flow so that the reparametrized flow is the magnetic 
flow ip of the constant curvature metric g\ and magnetic field given by a constant 
multiple of the area form of g\. Hence the flow ip is algebraic and thus it has smooth 
weak foliations. Consequently, <fi has also smooth weak foliations. 

Suppose now that <fi has both weak foliations of class C 1,1 . By f[JJ Theorem 4.6] we 
can perform a smooth time change on so that the new flow is smoothly conjugate to 
a quasi-fuchsian flow as defined in ^2]- Since the Godbillon-Vey invariant of the weak 
foliations of a quasi-fuchsian flow equals 4n 2 x{M) and the Godbillon-Vey invariant 
does not change under conjugacies we conclude that for our thermostat we must 
have 

gv ± = 4tt 2 X (M). 

By Corollary 13.41 this implies that 

V^ + 9) = 0. 

If we let h = r + 9, we can think of h as smooth functions defined on M. By 
Lemma \2. 71 we have 

9 = F(log(/i+ - h~)) + h + + hT. 

But if g is a smooth function on M, then F(g on) = dg and the last equality implies 
right away that h + + h~ = and 

9 = F(log(/i + - hr)). 

Thus 9 is an exact 1-form (i.e. e = —VU for some smooth function U) and by Lemma 
12. 21 we may perform a time change so that we just have to deal with a magnetic flow 
whose functions r ± satisfy Vir^ = 0. Again, the equality in Lemma \2. 71 tells us that 
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r + + r =0 and F(log(2r + )) = and thus r 1 * 1 are constant functions. If we now input 
this information into the Riccati type equation from Lemma [2.51 we obtain 

H(fon) = r 2 + K + f 

which implies that / and K are constant functions as desired. 

□ 

5. Proof of Theorem B 

We first characterize the thermostats for which the Godbillon-Vey invariant is max- 
imal. 

Proposition 5.1. Let M be a closed oriented surface and let <ft be an Anosov thermo- 
stat with divE = 0. Suppose the weak foliations are of class C 1,a with a > 1/2. Then 
there is a weak foliation with gv = Att 2 x{M) if and only if g has constant curvature 
—c 2 and there exists h 6 C°°(M, K) such that h 2 + f 2 = c 2 and 

e = \i(fVh-hVf). 

To prove this proposition we need some preparations. Following V. Guillemin and 
D. Kazhdan in [T3] we introduce the following first order differential operators: 

V+ :=(X-iH)/2 

and 

:= (X + iH)/2. 

Let L 2 (SM) be the space of square integrable functions with respect to the Liouville 
measure of SM. We now summarize some of the main properties of these operators 
(cf. d). 

(1) L 2 (SM) decomposes into an orthogonal direct sum of subspaces H n , nGZ, 
such that on H n , — i V is n times the identity operator; 

(2) r) + extends to a densely defined operator from H n to H n+ \ for all n. Moreover, 
its transpose is — r/_; 

(3) let C^(SM) = H n nC°°{SM). The operators ij± : C~ -> C™ ±1 are first order 
elliptic differential operators. 

Given a smooth 1-form 9 we can decompose 9 as 

9 = 0_! + 9 1 

where 

2fl_! =9 + iV(9), 
19 x =9-iV{9). 

Clearly 9 ±1 e H ±1 . Note that rj + 9^i = 

The following easy lemma will be useful later on. 
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Lemma 5.2. The form 9 is closed if and only if Qrj^Oi = 0. The form 9 is coclosed 
if and only ifdlr)-9i = 0. Also, 9 is closed if and only V(9) coclosed (and hence 9 is 
coclosed if and only ifV(9) is closed since V 2 {9) = —9). 

Proof. Note that 43fy_0 a = X{9) + HV{9) and that 43??_fli = H{9) - XV{9). Let e 
be the vector field dual to 9. Then x (v) = (e(x),v) and V{9) = (e(x),iv). We now 
compute 

X(9)(x,v) = (V v e,v). 
Using the expression for the flow of H given in Section 2 we also get 

H(9)(x,v) = (V lv e,v). 

Thus 

X{9) + HV{9) = (V„e, v) + (V„e, iv) = dive 

and 

H(9)-XV{9) = (V iv e,v) - (V v e,iv) = -d9(v,iv). 

□ 

Proof of Proposition [Ql Suppose first there is a weak foliation with gv = 4tt 2 x{M). 
By Corollary 13.41 there exists a function h G C 1 (M, M.) such that r = h o n — 9 is a 
solution of the Riccati type equation of Lemma 12.51 (in what follows we write h also 
for h o 7r and similarly for / and K) . Thus 

F(h) + h(h-9) + K- H(f + V{9)) + (f + V(9)) 2 - F(9) = 0. 

Using that F = X + (f + V(9))V and that V(h) = we obtain 

X{h) -h9- H{f) + fV{9) +K + h 2 + f - (X(9) + HV{9)) = 

and since dive = 0, 4^rj_9i = X(9) + HV(9) = and therefore 

(10) X{h) -h9- H{f) + fV{9) + K + h 2 + f 2 = 0. 

Now observe that X(h) - h9 - H(f) - fV{9) e E_ x © E x and K + h 2 + f 2 e H , so 
equation (JTUJ) is equivalent to the pair of equations 

(11) X(h) - h9 - H(f) + fV(9) = 0, 

(12) K + h 2 + f 2 = 0. 

Consider the complex valued function p := h + if and let uj := X(h) — h9 — H(f) + 
f'V(9). A simple calculation shows that 

V-(p) - 0-iP = (w + iV{u))/2 = 
It follows that equation (JTTJ) which is just u = 0, is equivalent to 

(13) 77_(p)-0_ip = O. 

Since ^± are elliptic operators, equation (fT^j) tells us that /i G C°°(M, R). 
Lemma 5.3. If p is a solution of M 6 J\) then 

V+V~{\P\ 2 ) = \V+(P) + 0iP\ 2 + 2\p\ 2 ^ V -0i- 
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Proof. This is an easy calculation in which one uses that r]_p = r] + p and that rj + r]_ = 
r]^r] + on H . We omit the details. □ 

Since we are assuming that dive = 0, Lemma 15.21 tells us that dtrj-Oi = and 
therefore, on account of the last lemma 

v+v-W) = \v+(p) + 0ip\ 2 >o. 

Since i] + r]_ on H is nothing but the Laplace-Beltrami operator in disguise, the last 
inequality implies that \p\ 2 must be a constant function on M and 

(14) V+ ( P ) + e lP = o. 

The fact that \p\ 2 is constant, combined with (|12|) shows that the metric g must have 
constant curvature, let us say, — c 2 . If we now combine equations (fT^j) and we 
obtain: 

X(p) - ipV(B) = 

which is equivalent to 

x(h) + fv(e) = o 

and 

X(f) - hV(B) = 0. 

Solving these equations for 9 we arrive at 

c 2 6 = fH(h) - hH(f) 

which shows that the conditions in the proposition are necessary to have a maximal 
Godbillon-Vey invariant. 

To show that the conditions are also sufficient we only need to observe that if the 
metric has constant curvature — c 2 and h is a smooth function such that h 2 + f 2 = c 2 
and 

c 2 6 = fH(h) - hH(f) 

then h and / satisfy the pair of equations (jllj) and (J12)) and thus if we let r = h — 6, 
then r is a C°° function which satisfies the Riccati type equation of Lemma 12.51 It 
follows that the smooth vector fields F and H + rV span a two dimensional bundle 
invariant under the Anosov thermostat. This bundle can only be one of the weak 
bundles which shows that there is a C°° weak foliation which has gv = 47r 2 x{M) by 
Corollary E31 

□ 

5.1. Proof of Theorem B. Let us show that (1) implies (3) and let J 7 be a weak 
foliation of class C 1 ' 1 . By [T3J Theorem 4.1], T is transversally projective and thus T 
is C 1 -conjugate to the weak foliation of a geodesic flow of constant negative curvature 
(see pages 178 and 179 in ^3])- As in the proof of Theorem A we can conclude that 
gv(T) = An 2 x{M) and Proposition (|5.1j) implies that (3) holds. On the other hand 
(3) implies (2): it suffices to check that the argument given at the end of the proof 
of Proposition (J5.1|) shows that if (3) holds then there is a weak bundle which is C°° . 
Since (2) obviously implies (1), the proof of Theorem B is now complete. 
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6. Proof of Theorem C and complements 

6.1. Proof of Theorem C. Suppose that <fi has both weak foliations of class C 1 ' 1 . 
We now argue as in the proof of Theorem A and note that the presence of q does not 
really affect the argument that shows that 9 is exact. By ^3 Theorem 4.6] we can 
perform a smooth time change on (f) so that the new flow is smoothly conjugate to a 
quasi-fuchsian flow as defined in |HZ] . Since the Godbillon-Vey invariant of the weak 
foliations of a quasi-fuchsian flow equals 4n 2 x{M) and the Godbillon-Vey invariant 
does not change under conjugacies we conclude that for our thermostat we must 
have 

gv ± = 4tt 2 X (M). 

By Corollary 13.41 this implies that 

V(r ± + 9-V(q)/2) = 0. 

If we let = r^ 1 + 9 — V(q)/2, we can think of h^ 1 as smooth functions defined on 
M. Note that V(X) = —9 + V(q) and hence by Lemma [2.71 we have 

9 = F(log(/i + - h~)) + h + + hT. 

The last equality implies right away that h + + h~ = and 

9 = F(log(/i + - h~)). 

Thus 9 is an exact 1-form and since we are assuming that it is coclosed, it must vanish 
identically. It follows that h must be constant functions. 

Let us now use that r = h + V(q)/2, where h = h is a constant, is a solution of 
the Riccati type equation of Lemma f2.5l with A = / + q. We obtain: 

F(h - V{q)/2) + {h+ V{q)/2){h - V{q)/2) + K - H(f + q) + (/ + qf + F{V{q)) = 

where F = X + (/ + q)V . Using that h is a constant and V 2 (q) = —4q we derive: 

(15) XV(q)/2 -H(f + q) + K + h 2 + f 2 -q 2 - [V(q)] 2 /A = 0. 

Now observe that XV(q)/2-H(f+q) G H^®Hi and K+h 2 +f 2 -q 2 -[V(q)} 2 /A e H , 
so equation (fl3|) is equivalent to the pair of equations 

(16) XV(q)/2-H(f + q) = 0, 

(17) K + h 2 + f - q 2 - [V(q)] 2 /A = 0. 
It is straightforward to check that 

det g (q + fg) = f 2 -q 2 -{V(q)] 2 /4 

and thus equation (JT7j) can be rewritten as 

K + h 2 + det g (q + fg) = 0. 

To understand equation (J16)) we will use the following easy lemma. 

Lemma 6.1. Let q be a symmetric 2-tensor which we also view as a function on SM . 
Then 8 g q\ SM = X(q) + HV(q)/2. 
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Proof. By definition of 5 g : 

S g q(x,v) = (V v q)(v,v) + (V iv q)(iv,v). 
Let 7 be the geodesic with initial conditions (x,v). Then 



X(q){x,v) = j t 



? 7 (t)('V(*)»'V(*)) = (V v q)(v,v). 

t=0 

Note that 

V(q)(x,v) = 2q x (iv,v) 
hence using the expression for the flow of H given in Section 2 we also get 

HV(q)(x,v) = 2(V io q)(iv,v). 

□ 

By applying V we see that equation (fl"6|) is equivalent to 

X(q-f) + HV(q)/2 = 
and thus the lemma implies that (|TT)j) is equivalent to 

s g (q - fg) = o. 

Summarizing, we have shown that ()16j) and ()17j) are equivalent to 

(18) S g (q - fg) = 0, 

(19) K g + h 2 + det g (q + fg) = 0. 

This shows that the equations listed in Theorem C are necessary for having weak 
foliations of class C 1 ' 1 . To show that they are sufficient we note that if we let = 
±h + V(q)/2, then is a C°° function which satisfies the Riccati type equation 
of Lemma 12.51 It follows that the smooth vector fields F and H + r ± V r span a 
two dimensional bundle invariant under the Anosov thermostat. This bundle can 
only be one of the weak bundles which shows that the weak foliations are C°° (with 
gv = Att 2 x{M) by Corollary IH.4jl . This finishes the proof of Theorem C. 

□ 

6.2. The space of solutions when / = 0. In this subsection we consider the 
solutions of the equations 

(20) 6 g (q) = 0, 

(21) K g + h 2 + det g (q) = 0, 

which arise from Theorem C when e = / = 0. 

Given a C°° metric g on M, its conformal class contains a unique metric go with con- 
stant curvature —1. The metric go determines a complex structure and the equation 
5 g (q) = simply says that q is the real part of a holomorphic quadratic differential. 
Write g = e 2u go and note that 

det g (q) = e~ 4u det go (q). 
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The curvature K g can be expressed in terms of u as 

K g = e- 2u (A go u-l), 

where A go is the Laplace-Beltrami operator with respect to the metric go (with the 
sign chosen so that it is non- negative) . Thus equation (J2T]) can be written as: 

(22) A go u = l-h 2 e 2u -e- 2u det go (q). 

We now explain why (|2"2"|) has a unique smooth solution u for each h and q fixed. 
Since h 2 > and det go (q) < (recall that q has trace zero) we can find constants 
■u_ < u + such that 

1 -h 2 e 2u ~ -e~ 2u - det go (q) > 0, 
1 - h 2 e 2u+ - e' 2u+ det go (q) < 0. 

It is well known that under these conditions (cf. 1201) the semi- linear elliptic 
equation (}2"2"j) admits a C°° solution u with w_ < u < u + which must be unique by a 
maximum principle argument. 

For each nonzero h, let Sh be the space of solutions of (f20|) and pTf . The discussion 
above shows that Sh is parametrized by pairs (go, q) where go is a metric of constant 
curvature —1 and q is the real part of a quadratic differential of the complex structure 
determined by go- Of course, the diffeomorphism group V of M and its identity 
component T>o act on Sh- 

Ultimately we are interested in the flow <fi up to smooth time changes and smooth 
conjugacy. Having this in mind, we note that we only need to consider S\. Indeed, 
if (g,q,h) is a solution of (j2"Uj) and (|2*T|) then (h 2 g,hq, 1) is also a solution and the 
corresponding thermostats are the same up to a constant time change (cf. Lemma 

lOIl- 
Given (g,q) G <Si, we let r^ 1 = ±1 + V(q)/2. The smooth 2-dimensional bundles 
spanned by F and H + r V give rise to two smooth transversal foliations J r± . The 
foliations are invariant under the thermostat </> defined by (g,q)- By the results of 
Ghys jT^j, each J r± is smoothly conjugate to the weak foliation of the geodesic flow 
of a metric g± of constant curvature —1. Let T denote the Teichmuller space of the 
surface M and note that each g± determines an element [g±\ G T. Thus the map 
Si 3 (g, q) I— > ([#+], [g~]) E T x T induces a map 

G : S x /Vo ^ T x T 

and since Si/T> can be naturally identified with T(T), the tangent bundle of T, we 
get a map 

G : T(T) -4TxT. 

We do not really know that much about G. It clearly maps the zero section of T(T) 
onto the diagonal of T x T and it is easy to see that if G([g ], q) = ([g+], [9-]), then 
G([<7o], -?) = ([<?-], [<?+])■ 
Question. Is G a diffeomorphism? 

In the next subsection we will see that G also takes values outside the diagonal of 
T x T. 
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6.3. Entropy production. Consider the thermostat (ft determined by an arbitrary 
function A £ C°°(SM). The next result is taken from [3] and we include a proof for 
completeness sake. 

Theorem 6.2. LetpEC°°(SM) be such that X(p) + HV(p)/k = for some positive 
integer k, and suppose 

K - H(X) + A 2 [(A; + l) 2 /{2k + 1)] < 0. 
Then there exists u £ C°°(SM) such that F(u) = p if and only if p = 0. 

Proof. We will use the following identity proved in |3J Equation (13)]. Given u £ 
C°°{SM) we have: 

(23) 2 / HuVFudn= / (Fu) 2 dfi + / {Hu) 2 dfi 

JSM JSM JSM 

- [ {K -H{\) + \ 2 ){Vu) 2 dfi. 

JSM 

Using that X{p) + HV{p)/k = and that H and X preserve the Liouville measure 
we obtain: 

/ HuV(p)dfi=— / uHV{p)dfi = k I uX(p)dfi = —k I X(u)pdfx. 

JSM JSM JSM JSM 

Since X{u) = p — XV (u) we derive 

HuVFudfi = —k I p 2 d[i + k XV(u)pdfi. 

SM JSM JSM 



Combining the last equality with (|23j) yields 

(2fc + l) / p 2 dfi-2k I XV{u)pdfi 

JSM JSM 

+ [ {Hufdn- [ (K - H(X) + A 2 )(VV) 2 d/i = 0. 

JSM JSM 

We may rewrite this equality as: 

f (vWTTp-^mYd, 

Jsm\ vwttj 



SM 



K-H W+ X><£± 



^)(V(u)) 2 dfi+ [ (Hu) 2 dfi = 0. 

1 / JSM 



Combining this equality with the hypotheses we obtain Hu = Vu = which implies 
right away that u must be constant. 

□ 

Remark 6.3. If p(x,v) = q x (v,...,v) where q is a symmetric fc-tensor, then the 
condition X(p) + HV(p)/k = is just saying that q has zero divergence (cf. proof of 
Lemma 16.1)1 . For such a p and k = 1 it suffices to assume that (ft is Anosov 3J . It is 
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unknown if the Anosov hypothesis is enough for k > 2. The problem is open even for 
geodesic flows. We refer to [21] for partial results in this direction when k = 2. 

We now note that if p G Fl-k © Hk, k > 2, then X(p) + HV(p)/k = is equivalent 
to saying that rj-Pk = 0, where p = p^k + Pk and 

p + iV(p)/k = 2p_ k , 
p - iV(p)/k = 2p k . 

But the kernel of the elliptic operator ^_ in C k ^'(SM) is a finite dimensional vector 
space which can be identified with the space of holomorphic sections of the k-th power 
of the canonical line bundle. By the Riemann-Roch theorem, this space has complex 
dimension (2k — l)(genus(M) — 1) (for k = 2 we get the holomorphic quadratic 
differentials). 

Suppose <p is Anosov and let p be the SRB measure. The entropy production of the 
state p is given by 

e^(p) := — J div F dp = — Lyapunov exponents 

where divF is the divergence of F with respect to any volume form in SM. Recall 
that V(A) is the divergence of F with respect to G. 

D. Ruelle has shown that e^(p) > and it is not hard to see (cf. PJ) that 
e <j>{p) — if and only if there is a smooth solution u to the coho mo logical equation 
F(u) = V(A). Equivalently <ft is volume preserving if and only if e^(p) = 0. Theorem 
16.21 gives a large class of thermostats with positive entropy production just by taking 
A = 9ft (pfc), where rj^p k = 0. In this case V(A) = $t(ikp k ). Note that is reversible 
if k is odd so we get lots of new examples to which the Fluctuation Theorem of G. 
Gallavotti and E.G.D. Cohen applies [3 00 CHI- 
Let (g, q) G S\ be a pair as in the previous subsection. Then V(q) is the real 
part of the restriction to SM of a holomorphic quadratic differential. If q is small 
enough then is an Anosov flow and by Theorem 16.21 is volume preserving if and 
only if q = 0. This shows that the map G above is interesting and that dissipative 
quasi-fuchsian flows do appear as thermostats associated with quadratic differentials 
as claimed in the Introduction. 
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